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Abstract 

Models of spontaneous wave function collapse modify the linear Schrodinger equation of standard 
Quantum Mechanics by adding stochastic non-linear terms to it. The aim of such models is to 
describe the quantum (linear) nature of microsystems along with the classical nature (violation 
of superposition principle) of macroscopic ones. The addition of such non-linear terms in the 
Schrodinger equation leads to non-conservation of energy of the system under consideration. Thus, 
a striking feature of collapse models is to heat non-relativistic particles with a constant rate. If such 
a process is physical, then it has the ability to perturb the well-understood thermal history of the 
universe. In this article we will try to investigate the impacts of such heating terms, according to 
the Continuous Spontaneous Localization (CSL) model, on standard evolution of non-relativistic 
matter and on the formation of CMBR. We will also put constraints on the CSL collapse rate 
A by considering that the standard evolution of non-relativistic matter is not hampered and the 
observed precise blackbody spectrum of CMBR would not get distorted (in the form of fi— type 
and y— type distortions) so as to violate the observed bounds. 
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1. INTRODUCTION 

Models of spontaneous wave function collapse jl, 2| aim to unify the dynamics of micro- 
scopic and macroscopic systems in order to answer the long prevailing question of Quan- 
tum "Measurement problem". The unification of microscopic (superposition of states) and 
macroscopic (violation of superposition principle) dynamics is accomplished by modifying 
the Schrodinger equation through adding non-linear stochastic terms. The non-linear terms 
in the modified Schrodinger equation ensures the breakdown of superposition principle at 
the macroscopic level and the stochastic nature of such dynamics indicates that the outcome 
of measurements would be probabilistic. The added non-linear terms act as amplification 
mechanism to ensure that these modifications have negligible impacts on microscopic system 
but are very efficient of localization for macroscopic ones. It is also important to note that 
the added stochastic terms also respect causality of the dynamics as deterministic non-linear 
evolution of Schrodinger equation leads to violation of relativity as shown in [3|. Among the 
many attempts of constructing such a modified dynamics, Quantum Mechanics of Sponta- 
neous Localization model (QMSL, later known as GRW model after the name of the authors 
4|), Quantum Mechanics with Universal Position Localization model (QMUPL 5], Q]) and 
Continuous Spontaneous Localization (CSL [7, 8]) model are worth noting. In some models 
it is considered that localization of wave-packets is a consequence of gravitational effects 



In this article we will concentrate only on CSL model and analyze the signatures of 



such a model while applying it in cosmology. 

According to the present understanding of the CSL model, a theoretical origin of the 
free parameters introduced in the scheme, i.e. the collapse rate A and the width of the 
localization r c , are yet to be determined. Hence values of such free parameters have to be 
obtained phenomenologically. Several proposed bounds on the free parameter A (considering 
r c w 10 -5 cm) can be found in j^, [l2 ]. 

One distinguishing feature of the CSL model, common to most collapse models, is that 
due to the presence of non-linear stochastic terms in the Schrodinger equation the total 
energy of the system does not remain constant and the non-relativistic massive particles 
within a system gain energy with a constant rate proportional to the CSL parameter A. 
Such heating effects are interesting to study for systems which are generally in thermal 
equilibrium. Study of atomic and nuclear systems [8] which demands the non-dissociation 
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of cosmic hydrogen by CSL heating obtains an upper bound as A < 1 s 1 , whereas studies 



of proton decay considering CSL models 



13] also lead to similar bounds. Comparing with 



the experimental bounds on photon emission from Germanium one can constrain the CSL 
parameter as A < 10~ n s _1 . Furthermore, in a cosmological scenario heating of protons 
through CSL mechanism over universe's life-time and the rise of CMBR temperature 

due to interactions of CMBR photons with these heated protons 13J have been studied to 
suggest a bound as A < 10~ 5 s _1 . Also thermal equilibrium of Inter Galactic Medium (IGM) 



131 ] has been studied to put an upper bound as A < 10~ 8 s~ 4 . 

However, as the thermal evolution of our universe is well studied and constrained observa- 
tionally, there can be other cosmological scenarios, apart from the evolution of IGM, where 
the imprints of the anomalous heating of CSL model can be studied. One such scenario, 
which we will exploit in the present article, is the thermodynamic equilibrium of matter 
and radiation before the formation of CMBR. Injection of energy during such epochs per- 
turbs the thermodynamic equilibrium of matter and radiation in the cosmic plasma which 
leads to distortions in the well-measured blackbody spectrum of CMBR. Such distortions 
are stringently constrained by observations. Here we will address two types of distortions 
of CMBR : the fi— distortion, which results from energy releases during the redshift span 
2 x 10 6 > z > 5 x 10 4 , and the y— distortion, which comes from energy releases during 
the redshift span 1100 < z < 5 x 10 4 . fi— distortion yields a non- vanishing, frequency- 
dependent small chemical potential of the cosmic photons which leads to a Bose-Einstein 
distribution in the high-frequency regime of the spectrum rather than the pure Planckian 
spectrum of a blackbody |l4j. Observations of COBE/FIRAS pointed out that such a dis- 
tortion should be very small and puts an upper bound on such distortions as \i < 9 x 10~ 5 



with 95% confidence 



151 ]. whereas an upcoming experiment PIXIE can probe such distor- 



tions up to \i ~ 5 x 10 [16] . On the other hand, the y— distortion of the spectrum is 
characterized by the lower density of photons in the low-frequency regime and increment 
of photon number in the high-frequencies with respect to a standard Planckian blackbody 
spectrum [ljj]. COBE/FIRAS observations put an upper bound on observed y— parameter 
as y < 1.5 x 10 -5 with 95% confidence jl5|. whereas PIXIE can put an upper bound on such 



distortion as y < 10 8 [16]. 



In this present article we will explore the possibilities of generating spectral distortions 
of CMBR due to CSL heating of non-relativistic particles. Before any such endeavor, one 
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should confirm that evolution of non-relativistic matter during the radiation and matter 
dominated era should not get affected by such anomalous heating of CSL and they evolve 
according to the standard cosmology. Then we note here that as such an anomalous CSL 
heating of non-relativistic matter can disturb the thermodynamic equilibrium of matter and 
radiation before CMBR formation, it can lead to spectral distortions of the CMBR where 
the amount of distortions will depend upon the strength of the heating of the particles. 
Thus by quantifying the amount of distortions such CSL heating can generate, one can put 
bounds on the free parameter A of the theory. 

However, in the present article we will not consider the effects of CSL heating on the 
evolution of relativistic matter including radiation. It has been a challenge to consistently 
formulate any relativistic generalization of the collapse models such as CSL and various 
attempts have been made in that direction. For a brief review of the attempts made to make 
a relativistic generalization of Spontaneous Collapse models one may refer to Thus, here 
we will assume that the relativistic matter evolve according to the standard cosmology. Any 
properly developed relativistic collapse model can, in principle, leave its imprint on CMBR 
and will be worth studying in future. 

According to standard Friedmann-Robertson- Walker (FRW) cosmology at Zi ~ 2 x 10 6 
the temperature of the cosmic soup is Tj 470 eV which is much lower than the rest mass 
energy of electrons which is 0.5 MeV. Thus the electrons and protons present in the cosmic 
soup during the concerned epochs i.e. 2 x 10 6 > z > 1100 can be treated as non-relativistic 
particles. Furthermore, if one tracks the density profile n(z) = no(l + z) 3 , where no is the 
present number density of baryons today, then one obtains the number of baryons per unit 
correlation cell of size rf. ~ 10 -21 m 3 (prescribed in CSL model) to be smaller than unity 
at Zi ~ 10 6 which falls off even further with the expansion of the universe. This indicates 
that considering the electrons and protons as free particles during such epochs is a good 
approximation. Hence the above discussion shows that during such epochs CSL model is 
applicable considering the massive particles in the cosmic soup to be free and non-relativistic. 
This enforces one to seek for possible spectral distortions in the CMBR generated due to 
CSL heating and to put bounds on A by quantifying the amount of distortions that can be 
yielded by such heating. 

We have organized the present article in the following manner. In Sec. (J2J) we will briefly 
review the CSL model and estimate the rate of heating of non-relativistic particles due 
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to presence of non-linear stochastic terms in the Schrodinger equation in the CSL sce- 
nario. Here, we will discuss two scenarios of CSL model : firstly the generic case where the 
strength parameter A is independent of mass and a variant scenario where A is dependent 
on the mass of the non-relativistic particles under consideration. We have mentioned be- 
fore that one should ensure that CSL heating would not disturb the standard evolution of 
the non-relativistic matter during radiation and matter dominated era. Sec. ([3]) is devoted 
to such discussions and to constrain the strength parameter A by demanding the standard 
cosmological evolution of non-relativistic matter throughout. In Sec. (jl]) we will first give 
a brief review of \i— distortion and then we will put bounds on both the scenarios of CSL 
model discussed above. Sec. §5§ is focused to discuss the y— distortion that CSL heating of 
non-relativistic particles can generate and constraining the CSL parameter from such dis- 
tortion of CMBR spectrum. In the penultimate section (Sec. (jSJ)) we will briefly consider 
the case where spontaneous photon emission is taken into account within the arena of CSL 
model and show that such feeble emission of photon in the high-frequency regime would be 
insufficient to generate any considerable distortion in the CMBR spectrum and thus would 
not be an important feature to constrain the strength parameter A. In the last section we 
will discuss all the bounds we will obtain throughout and then conclude. 



2. REVIEWING CSL MODEL AND EXPLAINING THE HEATING OF PARTI- 
CLES 

In this section we will briefly review the CSL mechanism and how it leads to heat- 
ing of non-relativistic particles with a constant rate. We will initially describe the GRW 
model of spontaneous collapse [4] where randomly acting non-linear terms are added to the 
Schrodinger equation, based on the assumption that the constituents of a system suddenly 
collapse to a localized wave-function with an appropriate rate. Concepts of such models are 
easy to conceive and the main result for energy gain by the non-relativistic particles are 
the same in both GRW and CSL models which makes the GRW model worth describing at 
this point. In literature, at later times GRW model has been extended to CSL model by 
appropriately introducing stochastic terms in the Schrodinger equation combined with the 
theory of continuous Markov processes in Hilbert space. We will discuss the CSL model in 
the later subsection. 
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2.1. GRW model of spontaneous wave function collapse 
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A deviant of quantum mechanics, called the GRW model |4j, allows for spontaneous 
localization of a particle with a mean rate A, where the evolution of the system between 
two successive localizations is according to the standard Schrodinger equation. According 
to this theory a localization operator, 
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-IS-M, (i) 

7T/ 

localizes a wave-function \ip) in space (at point x) yielding IV'x) II' "where |-?/> x ) = L^i/j) 

is the wave-function of the i th particle localized at x. Here 1/y/a = r c quantifies the accuracy 
of the localization at position qj. Due to the stochasticity of the spontaneous localization, 
a system in a pure state, which consists of a single particle, is transformed into a statistical 
mixture : 

|^>M-> J d 3 xLi\iP){^\Li = T[\i;){i;\}. (2) 

Analogously, if the system is initially in a statistical mixture, given by the operator p, 
then the localization of the system is determined by T[p] with a probability A. Hence the 
statistical state p after time-interval dt can be obtained as 



p(t + dt) = (1 - Xdt) 



p(t)-~[H,p(t)]dt 



XdtT[p(t)}, (3) 



where the first term on R.H.S. indicates that the system evolves according to Schrodinger 
equation if it is not localized within the interval dt. Thus the evolution of the system, 
represented by operator p, will be according to the master equation : 

j t p(t) = ~[H,p(t)]-\(p(t)-T[p(t)\). (4) 



A simple example of free particle Hamiltonian in one dimension has been considered in 
where the solution of the master equation in coordinate space has been obtained as 



4] 



{q ' lp{t)lq " ) = Mj dk J d y^ kVF (k,q'-q",t)( q ' + y\p s (t)\ q " + y} J (5) 
where (q' |ps(^)| q") is the solution of the pure Schrodinger equation and the factor 

F(q,k,t) = e- xt+x ti d ^ (q - kT/m)2 (6) 



6 



encapsulates all the dynamics of localization. A — > yields F(q, k, t) = 1 which shows in this 
limit (q' |ps(^)| q") = (q' \p{t) \ q"), as expected. The dynamical evolution of this free particle 
system yields a spread in the momentum as [4| 

(f) = tr[p 2 p(t)] = (P 2 )s + ^-t, (7) 
which can be derived using Eq. (j3J). Here, {p 2 )s is the conserved momentum for the free 

2 

Schrodinger evolution. Thus for a non-relativistic particle with energy E = 2— the spread 
in energy will be 

{E) ^E )s + ^ t , (8) 

which shows non-conservation of energy of the system with time. Therefore, one can infer 
that due to spontaneous localization of particles in a system a non-relativistic particle gains 
energy with a rate 

dE_ _ 3a\h 2 

~dt ' ~ Am ^ ' 

This heating rate has been obtained for one non-relativistic particle. To investigate effects 
of such spontaneous localizations in cosmology one has to obtain the total heating rate for 
the cosmic plasma containing non-relativistic electrons and protons. 



2.2. Continuous Spontaneous Localization (CSL) model 



Out of many spontaneous collapse models, the one that is commonly used in physical 
applications is the Continuous Spontaneous Localization (CSL) model which generalizes the 
original GRW model 4j] to systems of identical particles. In the CSL model the collapse 
of the wave-function happens continuously in time in contrast to the GRW model where 
the collapse of wave-function happens discretely. The parameters introduced in this model, 
in a similar way as had been done in GRW model, are the coupling constant 7 and the 
correlation length r c (which is conventionally taken as 10 -5 cm). The modified Schrodinger 
equation takes the form in CSL model lj 



-Hdt + y/j I rf 3 x(iV(x) - (JV(x)) t )dWi(x) 



-1 J rf 3 x(iV(x) - (N(x)) t ) 2 dt 



IV*>- 



(10) 



Here the operator iV(x) is an averaged density operator given as 

iV(x) = I d 3 yg(y - *y (y, s)a(y, s), (11) 

s J 

where the sum is over various particle species s, with a|(y)a s (y) being the number density 
operator of species s and W t (x) is the family of standard Wiener process for each point 
in space. Moreover, g(x) is a spherically symmetric, positive, real function peaked around 
x = with the normalization J d 3 xg(x.) = 1 and can be written as 

9M = ^fr^ lTl - (12) 

The collapse rate A of GRW model is related to that of the CSL parameter 7 as 



A "8^7V^ (13) 



The corresponding master equation in CSL model as for Eq. (j3j) in GRW is given by [l|: 



dt h 



[H,p(t)]+ 1 J d 3 xiV(x)p(t)iV(x) - I J d 3 x{N 2 (x),p(t)}. (14) 



It is to be noted from this relation that A of GRW model is directly proportional to the 
parameter 7 of CSL model apart from some constant factors. In such bound on A 

will directly indicate a bound on 7. Hence we will consider bounds in A for further discussions 
and those can be converted into bounds on 7. 

In a generic CSL model scenario the introduced parameter 7 (or, equivalently, A) is 
independent of the mass of the constituent particles. In literature, a variant of such a model 
is also discussed where the parameter A is mass- dependent as 

A(m) = A f— ) , (15) 
\m N J 

where is the mass of a nucleon. Such a model is motivated from the feature that in 



,8 



this scenario the collapse rate is different for different species bearing different masses 
We provide a justification for such a choice. The calculations of the reduction rate of wave- 
functions in |l[ shows that the off-diagonal elements of coordinate space density matrix for 
a single nucleon approach zero exponentially with a reduction rate given by 

T R = A (l - e^ 2 /^) , (16) 
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which for |x| > r c becomes 

T R w A. (17) 

For n nucleons within a radius smaller than the correlation length, in a superposition of 
states with distance larger than r c , this rate is multiplied by n 2 to yield 



R 



n 2 \. (18) 



In mass-proportional CSL model one wants the collapse rate of a massive single particle 
composed of n fundamental particles to be the same as the collapse rate of n single funda- 
mental particles. This is true when the strength parameter A is quadratically proportional 
to the mass of the constituent particle i.e. (3 = 2. In that case the reduction rate of the 
system of n fundamental particles becomes 

T R = Xon 2 (—) 2 . (19) 
\m N J 

Thus in this article, along with the mass-independent CSL model with strength parameter 
A, we will also analyze the mass-dependent CSL model where the strength parameter Ao is 
for the quadratically mass-dependent case and compare the bounds obtained in both the 
cases. 

n 

It has been shown in [8] that due to the presence of non-linear terms in the Schrodinger 
equation in the CSL model the energy will not be conserved (as has also been discussed for 
GRW model). Since, as discussed in the introduction, the number of baryons in the universe 
per unit correlation cell of size r 3 ~ 10~ 21 m 3 is smaller than unity, for the considered times, 
the effect of the identity of particles can be neglected and the rate energy increase in such a 
scenario is similar to that in the GRW model which is given in Eq. (Q. For CSL model the 
non-conservation of energy is quantified as jg] 

(E) = 3 -^t. (20) 
4m 

Such a heating rate has been obtained for one non-relativistic particle. To investigate effects 
of such spontaneous localization in cosmology one has to obtain the total heating rate of 
the cosmic plasma containing non-relativistic electrons and protons. Given the number 
density n s for each species (s = e,p) the total energy density gain by each species can be 
straightforwardly obtained from the above equation as 

cte, 3a\h 2 



dt 4m, 



-n s , (21) 



taking into account that during the epochs one is interested in, these particles are non- 
relativistic and behave as free particles in the cosmic plasma. 

While applying a generic CSL model in cosmology, where the parameter introduced in 
the theory is mass- independent, the contribution to the total change in energy density will 
come from the electron fluid as — ^> — and n e ~ n„. Thus in such a scenario the change 

m e m p e P ° 

in total energy can be written using the above equation as 

de 3a\h 2 , , 

«=*sr n - (22) 

On the other hand, if one considers the variant scenario of CSL model where the parameter 
A is dependent on mass of the particle quadratically, then the proton fluid will contribute 
more to the total change in the energy density and in such a case one can write 

de 3aAoft 2 



dt Am p 



-n p , (23) 



where the mass of a nucleon can be taken as that of a proton m p . 

It is worth pointing out at this point that the origin of the stochastic field Wt(x) in 
the non-linear Schrodinger equation of CSL model is not yet known and thus it is difficult 
to point out at this moment how such a stochastic field can be cosmologically accounted 
for. In this work, we will thus not consider this stochastic field at all but only analyze the 
energy-increase of standard non-relativistic particles while interacting with such fields. 



3. BOUNDS ON CSL PARAMETER FROM RADIATION DOMINATED AND 
MATTER DOMINATED ERA 

Before investigating the spectral distortions in CMBR caused by CSL heating of non- 
relativistic particles, it is important to ensure that such heating does not hamper the stan- 
dard evolution history of the universe. We can only consider those epochs where the particles 
can be treated as non-relativistic particles and hence allow one to apply the methods of CSL 
model. 

In standard FRW cosmology the scale factor a is related to the corresponding redshift z 
as 1 + z = ^ where a is the present scale factor. This implies that ^ = — It is also 
considered that the universe evolves adiabatically which implies that the comoving entropy 
is always conserved. For such isentropic processes a 3 T 3 remains constant in a comoving 
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volume which yields T oc - oc (1 + z) i.e. T(z) = T (l + z) where T = 2.73 K is the present 
temperature of the CMBR. Writing the temperature in eV units one has 

T(z) « 2.4 x 10~ 4 (l + z)eV, (24) 

which implies that when the temperature of the universe is 0.5 MeV (i.e. of the order of 
the rest mass energy of the electrons) the corresponding redshift would be 2 x 10 9 . Thus for 
radiation dominated (RD) universe we would be interested in the redshift span of 2 x 10 9 < 
z < 3233 and for that of matter dominated (MD) era would be 3233 < z < 1 (where 
z ~ 3233 is the epoch of matter-radiation equality). 

Due to CSL heating the energy density of matter will evolve differently than standard 
cosmology. According to standard cosmology the energy density of matter gets diluted with 
the expansion of the universe as 

^ = -3H(z) PM (z), (25) 

where H = - is the Hubble parameter and we have taken the pressure p = for matter 
(non-relativistic particles). But the change in the matter energy density with time (^f-) 
in this case is affected by both the expansion of the universe and the CSL heating which 
is proportional to the matter number density present during that time in the cosmic soup. 
Then one can write the evolution of the matter density using Eq. ( 121 p . in the case where 
radiation and matter are not tightly coupled, as 

dt = -3H{z)pm{z) + 2^ n M {z). (26) 

s=e,p s 



In standard cosmology the energy density pu of non-relativistic electrons and protons can 
be written in terms of their number density nyi as 

Pm = n M (m e + m p )c 2 , (27) 

where the number density of electron n e and that of proton n p are the same n e ^ n v ^ um- 
Using this in Eq. ( 12 6 p one can write as 

^# = -3^ + Y 3Xh2a , 9 dt. (28) 
Pm{z) a ^ Am s (m e + m p )c 2 

With the definition of Hubble parameter H and the relation 1 + z — — it can be seen that 

dz . „ N 

* = (29) 
11 



and along with the relation 
written as 

dp M (z) 



da 
a 



the evolution equation of energy density can be 



dz 



Pm{z) 



E 



s=e,p 



3\h 2 a 



dz 



Am s (m e + m p )c 2 (1 + z)H(z) 



(30) 



The Hubble parameter H(z) = - can be written in terms of dimensionless quantities using 
the Friedmann equation as 



A ) 



(31) 



where f2j = pi/(3rripH 2 ) is called the density parameter for the matter species i and Qr , 
^m ; ^fc an d ^a are the present radiation, matter, curvature and dark energy density 
parameter respectively. During RD era, the radiation fluid contributes dominantly to the 
total density and thus the Hubble parameter can be written as 



H(z) = H Ql o (l + z 



(32) 



whereas during MD era that would be 



H(z) 



(33) 



Let us consider the RD era first. It is worth noting at this point that during radiation 
era matter and radiation will be tightly coupled and during this period Compton scattering 
is much more efficient than CSL heating (as will be argued below). Also as the ratio of 
entropy in baryons and photons is 10~ 9 the baryons will lose most of its energy, gained by 
CSL heating, to photons via Compton scattering retaining only a fraction / of the order 
10~ 9 . Keeping this in mind and using Eq. fl32|) in Eq. fl30|) one gets 



dp M {z) 



dz 
'l + z 



Pm(z) 

solving which one gets 

p M} = exp <J i^CSL.RD 

where we have defined 



-E 



3f\h 2 a 



dz 



(34) 



i + ^ 



l + Zi 



PMt 



CSL,RD 



E 



3 x 10- 9 \h 2 a 



s =e,p 8m s (m e + m p )c 2 VL 2 RQ HQ 
12 



(35) 



(36) 



It is evident from the above equation that in absence of CSL heating (i.e. with A = 0) we 
have the standard cosmological evolution of matter density as 

07) 

It is important to note that before recombination the free electrons and protons are the 
dominant component which contribute to the total energy density of the universe as non- 
relativistic particles i.e. p m ~ p e + p P - As the number density of electron and protons are the 
same, as has been considered earlier, and the protons are much heavier than the electrons, 
one can consider that the non-relativistic fluid of the cosmic plasma is dominated by the 
proton fluid (p^/ ~ Pp)- Following such arguments one can ignore the electron contributions 
in Eq. ( |35l) and write 

3x l(T 9 Aft 2 a 

^CSL.RD ~ Yl 1 • ( 38 ) 

s =e, P 8m s m p c 2 {l'x o H 

Thus to ensure that the standard evolution of the non-relativistic plasma during RD era is 
not disturbed due to CSL heating one requires 

1 1 



-KcSL.RD 



< 1, (39) 



where Z{ « 2 x 10 9 and Zf « 3233. 

Thus in the generic scenario of CSL model (where the total change in energy density will 
depend upon the change in the electrons' energy density) one obtains the constraint on A 
using Eq. ([2"2]) as 

A| RD < 5 x 10 10 s-\ (40) 

while for the variant scenario with A quadratically proportional to the mass one puts an 
upper bound on the parameter using Eq. ( 123"]) as 

A | RD «10 14 s- 1 . (41) 

During MD era using Eq. fl3"3"|) in Eq. (150]) one gets 

dp M {z) _ 3 dz _ f\h 2 a dz ^ 

Pm( z ) 1 + z s=e>p 4ms ( me + mp )c 2 H nl Io (1 + ' 
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and solving which yields 

p Mf = exp | i^CSL,MD 
where we have defined 

-Kcsl.md = 



1 + Zf)* (1 + Z t ) 




(43) 



E 



f\H 2 a 



(44) 



s=e, P Qm s (m e + m p )c 2 VL 2 MQ HQ 

But the situation changes in the MD era i.e. just before and after the recombination. 
Here we will consider the redshift span of 3233 < z < 1. The recombination happens at 
z ~ 1100 which implies that for z < 1100 the matter and radiation will seize to remain 
tightly coupled and thereafter they will hardly interact. Thus, nearly all the heat gained 
by the non-relativistic matter due to CSL heating after z ~ 1100 will be retained in the 
matter itself, yielding / ps 1. Thus the strongest bound on A will come from the era 
1100 < z < 1. Also after recombination the energy density would be dominated by neutral 
hydrogen (pm ~ Ph)- Taking pu ps p v we can consider pu ~ p P throughout 1100 < z < 1. 
Following such arguments then one can drop the contribution of the electrons in the MD 
era as well and defining 

v - \h 2 a 



CSL,MD ~ 1 ' 

s =e,p Qm s m p c 2 Vt 2 MQ H 



(45) 



one should have 



CSL,MD 



«1, 



(46) 



[1 + Zf)* (1 + ZiY- 

to ensure that the non-relativistic fluid evolves according to the standard cosmology during 
the MD era. Taking Zi ~ 1100 and Zf ps 1 for the MD era the above condition puts an upper 
bound on the CSL parameter A for the generic case as 



-^Imd 



«4x 10"V 



(47) 



where we have used Mo ps 0.04 as the present baryon density parameter. Similarly for the 
mass-dependent case one gets the upper bound as 



Ao| M D<°- 7s l - 



(48) 
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4. BOUNDS ON A FROM ii-TYPE DISTORTION OF CMB 



The thermalization history of the CMBR photons is affected by any unusual energy in- 
jection during its period of thermalization with the baryons in cosmic soup and results in 
spectral distortions of the observed CMBR blackbody spectrum [jj], fljl]. The unusual energy 
release in these early e poc hs can be due to decay of relic particles [181 ]. by evaporation of 



primordial black holes [19j . WIMP annihilation 20], damping of sound waves in primordial 
plasma [2l| or by other astrophysical mechanisms as has been mentioned in [20]. Ther- 
malization of the photons requires non-conservation of photon number in the cosmic soup. 
Cosmologically double-Compton scattering and bremsstrahlung processes are the ones which 
produce photons and thus are very significant in the thermalization of photons in the cos- 
mic plasma. These photon number non-conserving processes are efficient in erasing possible 
spectral distortions in the early epochs. On the other hand, the photon-number conserving 
process of elastic Compton scattering, 7 + e~ — > 7 + e~, is the most dominant process which 
couples electrons with photons at early times and redistributes the photons in frequency. 
As the photon number non-conserving processes, mainly the double Compton scattering 
process, become inefficient by z ~ 10 6 and photon-number conserving Compton scattering 
process then dominates the thermalization process, any energy release after z < 10 6 leads 
to spectral distortion of CMBR. 

Energy releases during 2 x 10 6 > z > 5 x 10 4 lead to a kind of CMBR spectral distor- 
tion, called /i— type spectral distortion, which yields a non-vanishing frequency dependent 
chemical potential of the CMBR photons. This epoch 2 x 10 6 >z>5x 10 4 is known as the 
/i— era in the literature. Due to such energy releases the Planck spectrum of the thermalized 
)hoton density relaxes to a Bose-Einstein distribution with a non-zero chemical potential fi 



14J as 



Wpl = e e/«sT _ i — > U ^ = e (^)/l B T _ I ( 49 ) 

Absence of significant distortion in the CMBR spectrum from the Black Body distribution 



observed by COBE/FIRAS puts an upper limit on such distortion as // < 9 x 10 5 [15]. 



An upcoming space-based mission, called PIXIE, will constrain the /i— type distortion up to 



H ~ 5 x 10~ 8 [16]. 



An interesting feature of these processes is that while double-Compton scattering and 
bremsstrahlung processes become efficient as frequency decreases, elastic Compton scat- 
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tering process, on the other hand, is independent of the frequency of the photons. Thus 
at low redshifts when the Compton scattering is the dominant process of thermalization 
of the cosmic plasma, the photon number non-conserving double-Compton scattering and 
bremsstrahlung processes are still efficient at low frequencies and are capable of retaining 
the black body spectrum 22]. Also, the photons, which are generating at low frequencies 
and getting scattered to high frequency regime, are not sufficient to maintain a Planck spec- 
trum at high frequencies. Thus the elastic Compton scattering establishes a Bose-Einstein 
spectrum only at the high frequencies of the spectrum. 

The total //— distortion generated due to energy injections and injection of high-energy 



photons during 2 x 10 >2>5x 10 can be written as as 
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^ = -i_f 3 *E-4*ty (50) 



2.143 \ e n J 

where 5e and 6n are the total changes in energy and number density respectively with respect 
to a Planckian distribution. A brief derivation of the above formula is given in Appendix (IA"|) . 
This remarkable result shows that the total energy distortion of /z— type is independent of 
the actual form of the energy injection. It can also be seen from the above equation that 
direct heating of the electrons contributes in the same way as injection of high-frequency 
photons as has been pointed out in 22]. 

It is important to note here that Compton scattering during these epochs is very efficient 
to maintain the equilibrium between the electron and photon temperature and keeps the 
cosmic fluid tightly coupled. Thus the CSL heating rate should be much smaller than the 
Compton scattering rate such that the thermal equilibrium of the cosmic plasma should be 
retained. During these epochs equilibrium is attained within a time scale given as 22] 

tc = - 7.63 x 10 19 (1 + zY A s, (51) 

where e 7 is the energy density of photons at that redshift. Also, in the above equation the 
scattering of photons with electrons is only considered as that with protons is 10 6 times 
slower and e~ + 7 — > e~ + 7 would be the dominant process to keep the plasma in thermal 
equilibrium. Thus one would expect the CSL heating rate to be much smaller than the 
rate of attaining equilibrium (t^, 1 ) via Compton scattering till the recombination ends at 
z ~ 1100, which using Eq. ( 150]) weakly puts an upper bound on the CSL heating rate for 
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the generic case as 

^2xl0 3 s~\ (52) 



A|cs < 3a/i 2 c 



2~1100 

where we have used the values m e w 9.1 x 10~ 31 kg, a ~ 10 14 m -2 , c = 3 x 10 8 m s _1 and 
/i = 6.62 x 1CT 34 kg m 2 s _1 . Similar bounds in the variant scenario where A is quadratically 
mass-proportional would be 



Am p m e c 2 x 



In a generic scenario of CSL model 



wSxlO's" 1 . (53) 

z~1100 



, |4[ , non-relativistic particles gain energy with a rate 
A without changing the photon number density at higher frequencies. Here we will assume 
that the amount of energy gained by the electrons or protons due to CSL heating will be 
transferred to the photons completely as during these epochs the number density of photons 
exceeds that of the electrons by many orders of magnitude. Thus the total \i— distortion 
due to CSL heating will be 

x 3 fecSL 

^ = 2l43— (54) 
where we have neglected the last term in Eq. (|50|) as there is no significant photon production 
in a generic CSL scenario at high-frequencies. During the /i— era the total amount of energy 
gained by the non-relativistic fluid in the cosmic plasma can be determined by using Eq. (1211) . 
Using Eq. ( )29|) in Eq. (12 ip yields the rate of energy gain by the non-relativistic particles 
with redshift as 

8e _ v 3«A£ 2 (1 + z) 2 

s=e,p v ' 

where we have used the relation between the present number density of electrons n eo and 

that in a particular redshift n e (z) as n e (z) = n eo (l + z) 3 . During the radiation era, i.e. 

i 

before recombination, the Hubble parameter H(z) can be written as H(z) ~ fijLf/o(l + z ) 2 

where Qr ~ 6.5 x 10 -5 and Hq 2 x 10 -18 s^ 1 are the present radiation density parameter 

and Hubble parameter respectively. Putting this value for H(z) in the above equation yields 

_ _ 3a\h 2 n eo 
5z ~ 2s 4 m i ■ 

s=e,p "Ro-^O 

The unperturbed Planckian photon energy at z would be 

e(z) = ^(l + z)\ (57) 



17 



which will give the fractional change in the energy of a pure Planckian spectrum at redshift 

z as 

-?=-£ I"* (1 + ^ (58) 

s=e, P * \i Ro n aBi 

where T ~ 2.73 K is the present temperature of the CMBR. Integrating the above equation 
over Zi = 2 x 10 6 to zj = 5 x 10 4 one gets the total fractional energy gained by the electron 
or proton plasma due to CSL heating during the /i— era as 



E 



2 



e ^ lQm s 2 fj _ 7^4 

s=e,p Rq Q B 



1 1 



(59) 



;i + ^)3 (l + ^)3_ 

In a generic case, where the CSL parameter A is independent of mass, we have shown 
that the gain in the electron energy density will dominate and this will also be the total 
energy gained by the photons via Compton scattering through electrons. Hence the total 
/i— distortion generated due to CSL heating will be 



H ~ x 



aXh n eo c 



2.113 l6m enloHQaBT 4 



1 



~1.2xHT 6 As, (60) 



L(i + ^) 3 (i + Zl f 

where apart from the previously mentioned values of the parameters we have used n, 
0.246 m -3 . Using the COBE/FIRAS measurement of \i— distortion in the CMBR spectrum 



e ~ 



15[ one then can put an upper bound on the CSL heating parameter as 

A| m ,cobe/firas < 70s -1 , (61) 



whereas upcoming experiment PIXIE [16] can put a more stringent constraint on the CSL 
heating parameter as 

AUpixie<4x IO^s- 1 . (62) 

Similarly for the case where A is quadratically mass-proportional, photons will gain en- 
ergy from both the electrons and the protons via Compton scattering. Hence, the total 
/i— distortion can be calculated using Eq. (123"]) and Eq. ( 1571) as 



3 aXoh 2 
H ~ x x 



2.1 13 16m p 2 h^_t4 



6.8xlO~ iU As. (63) 



Thus bound on Ao from observations of COBE/FIRAS would be 

AoU,COBE/FIRAS ^ 10 5 S^ 1 , (64) 

and that of from PIXIE would be 

AoU,pixie<74s- 1 . (65) 
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5. BOUNDS ON A FROM y-TYPE DISTORTION OF CMB 



Energy release in the early universe during 2<5x 10 4 leads to a y— type CMBR spectral 



distortion 17j. At these low redshifts Compton scattering is not effective enough to maintain 
the full kinetic equilibrium between photons and electrons. Due to energy releases if the 
temperature of the electrons becomes greater than that of the photons then the low-energy 
photons will be upscattered in frequency via Compton scattering. Thus a y— type distortion 
is characterized by a deficit of photons in the low frequency regime with an increment 
of photons at high frequencies in comparison with a standard blackbody spectrum. The 
efficiency of Compton scattering process to maintain a blackbody spectrum is quantified by 
a parameter, called the Compton y— parameter, as 



V = I Kb(Te 2 Tl) N e a T cdt, (66) 
m e cr 



where Kb is the Boltzmann constant, T e and T 7 are the temperatures of the electrons and 
photons in the cosmic plasma respectively and N e is the electron number density. The present 
constraint on such distortion of CMBR comes from COBE/FIRAS experiments which puts 
an upper bound on ,-pa— a s y < 1.5 X 10- with 95% confidence Q. The upcouung 
experiment PIXIE can probe CMBR distortions up to y < 10 8 [161 ]. 

Energy injection of fe 7 in the photon energy density during z < 5 x 10 4 can yield a 



Compton y— parameter as 
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23] 



fy = if. (67) 



A brief derivation of the above formula is given in Appendix (JB]). In this work we are 
mainly interested in y— type distortions arising during the pre-recombinational era i.e. dur- 
ing 1100 < z < 5 x 10 4 . This is because any energy injection occurring after the completion 
of hydrogen recombination would not imprint any additional traces in the CMBR distor- 
tions. Later on, in post-recombinational epochs (z < 800), different physical mechanisms, 
like interaction of CMBR with hot intergalactic gas known as thermal Sunyaev-Zeldovich 



effect 



241 ]; supernova remnants at high redshifts or large-scale structure formation giving 



rise to shock waves 



2a 



261 ]. can also contribute to the y— distortion observed in CMBR. 
Contribution to y— distortion coming from such effects have been ignored in this study. 

As the dominant process during these epochs is still the elastic Compton scattering, 
we follow the same procedure to calculate the total energy gain by the photons, due to 
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CSL heating of non-relativistic electrons or protons present in the cosmic plasma, as we 
have discussed while calculating the \x— distortion. It is important to note here that during 
the y— era the redshift span 5 x 10 4 < z < 3233 is radiation dominated while the rest 
3233 < z < 1100 is matter dominated. Thus generalizing Eq. ([5T?j) and Eq. (|57p to reflect 
this fact and using those in Eq. (157)) in the redshift span Zi = 5 x 10 4 to Zf = 1100 one 
computes the total y— distortion generated due to generic CSL heating as 



1 3a\h 2 
~ 4 X 4m e 4a B T*H 



+ 



'l + zff 



1 + z, 



R\3 



(1 + z 



M ^ 



0.2 A s. 



(68) 



Here we have taken zf- w 5 x 10 4 and zf ~ 3233 for the radiation dominated epoch in 



y— era and similarly zf 1 ~ 3233 and z¥ 1100 for the matter dominated epoch. Thus 
the constraint on A from the measurement of y— parameter of CMBR by COBE/FIRAS 
experiment [15] would be 



A 



y,COBE/FIRAS 



< 8 x 10- & s- 1 , 



(69) 



whereas in future measurement obtained by PIXIE [16j can put an upper bound on A as 



A 



y.PIXIE 



< 5 x 10~ 8 s- 1 . 



(70) 



Similarly for the variant model where the parameter A is mass- dependent one has for the 
total y— distortion as 



y 



1 3aAo^ 2 n eo c 



x 



4 " 4m p 4a B T*H 



5 

+ — r 



1 



1 



[l + zf)i (l + zf) f 



;i + zf)* (i + zff 



~ 10~ 4 As. 



(71) 



Thus the upper bound on Ao obtained from COBE/FIRAS observations would be 

Ao|j/,COBE/FIRAS ^ 0.14 s" 1 , (72) 
and that coming from PIXIE's upcoming observations 



Ao 



j/,PIXIE 



< lO" 4 s 



4 „-l 



(73) 
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6. BRIEF DISCUSSION ON CSL MODEL WITH SPONTANEOUS PHOTON 
EMISSION AND ITS COSMOLOGICAL IMPLICATIONS 



In the literature it has been shown that for collapse models, atomic systems 13|, |27J, |28] 



or free electrons 



29] emit radiations spontaneously when the noise term in the evolution is 



treated perturbatively. The rate of creation of photons with wavenumber k due to sponta- 
neous collapse of a (free) particle species s can be calculated as 

dT(k) Xah 2 e 2 



dk 27T 2 eoc 3 m 2 k 

The change in the number density n 7 of photons due to this process would be 



dn 7 



Xah 2 e 2 n s0 (l + z) 3 log k 



1(T 48 log(fc)A(l + 2) 



(74) 



(75) 



csl 2n 2 e c 3 m 2 
where s has been taken as electrons. This injection of photons into the cosmic plasma 
will also contribute to the spectral distortion of CMBR in principle, as can be seen from 
Eq. ( 1501 . However, we can see from the above equation that the injection of photons due to 
spontaneous collapse will not provide any significant contribution at higher frequencies and 
will be suppressed by n as well. Even at lower frequencies this process is clearly subdominant 
compared to other photon producing processes like bremsstrahlung and double Compton 



scattering where the rates of photon creation are 



22|] 



~dt 



KT 18 log (2.25/x e ) 



and 



br 



(Iriry 



k 3 



[1 



l)]z 



5/2 



10 



-32 



DC 



A- 3 



-I(t)[l-n J (e B ' - l)]z 5 



respectively, for small x e . In the above equations 

hck 



k B Te 



with T e being the electron temperature and 

I(t) — I dx e x\ n 7 (n 7 + 1). 



(76) 



(77) 



(78) 



(79) 



We can see from above that for small x e and large n 1 



dt 



CSL 



10 30 A; 3 log(k)e x * 



dn-y 
~dT 



log(2.25/x e )[l -n 7 (e*« -1)]' 



,1/2 



(80) 



J br 
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and also 





Ally 






dt 


CSL 




dn^ 






dt 


DC 



Thus such a process at low frequencies are subdominant to other photon creating processes 
like double Compton scattering and bremsstrahlung which provide sufficient number of low 
energy photons to establish blackbody spectrum at lower frequencies. However, these photon 
creating processes including the one of spontaneous collapse model become more and more 
inefficient as frequency of the generated photons increases. Also the strength parameter A 
of the collapse model is stringently constrained to very low values (as low as 10~ 5 ) from 
present observations of CMBR spectral distortions. These arguments along with the fact 
that the spectral distortion is proportional to the number of photons generated which is 
suppressed by the total number of photons (the factor — of Eq. (150]) ) indicate that the 
generated photons at high-frequencies will not be able to yield significant distortions in the 
CMBR spectrum. Also — in the above equation can be significant if one takes A to be 
large which then violates the perturbative analysis of the model. So such a photon emitting 
collapse process will hardly improve the bounds on the strength parameter A obtained from 
ji— type and y— type distortions. 



7. DISCUSSION AND CONCLUSION 



In the present article we have addressed the effects of CSL heating of non-relativistic 
particles on standard cosmology and formation of CMBR. As the thermal evolution of our 
universe has been studied and understood rigorously by both theoretical and experimental 
means during the past half a century, any anomalous heating, which can disturb the thermal 
evolution of our universe, thus can be constrained stringently by observations. Motivated 
by this, we endeavor to investigate the effects a model like CSL, which heats the non- 
relativistic particles with a constant rate, can have on cosmological evolution and on CMBR 
formation. In the following Table U we summarize the bounds we obtain by ensuring that the 
CSL heating would not hamper the standard cosmological evolution of the non-relativistic 
particles and also it would not generate substantial distortions (fi— type and y— type) in 
the observed precise blackbody spectrum of CMBR. It is to be noted from Table |T] that 
the strongest bounds are coming from the observations of y— distortion in the CMBR. Such 
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TABLE I: Bounds on CSL strength parameter 



Case 


A (in s- 1 ) 


Ao (in s 1 ) 


Bounds from RD era of standard cosmology 


«5x 10 10 


< 10 14 


Bounds from MD era of standard cosmology 


«4x 1(T 4 


<C 0.7 


Bounds from comparing rates of Compton scattering and CSL heating 


2xl0 3 


3 x 10 6 


Bounds from COBE/FIRAS observation of ^—distortion 


70 


10 5 


Bounds from PIXIE future observation of fj,— distortion 


4 x 10~ 2 


74 


Bounds from COBE/FIRAS observation of y— distortion 


8 x 10~ 5 


0.14 


Bounds from PIXIE future observation of y— distortion 


5 x 10~ 8 


10" 4 



a bound on the strength parameter A or Ao is within the detection range of the proposed 
diffraction experiments from fullerene and larger molecules j^, [l^] . It suggests that apart 
from laboratory experiments, cosmological scenarios are also very important to test any 
consistent collapse model which is in accordance with the standard FRW framework of 



cosmo 
In 



ogy. 



131 ] a stronger bound on the strength parameter Aq has been obtained from the ob- 



served thermal equilibrium of IGM in the redshift span z ~ 6 — 2, by comparing the rate 
of heating of IGM through CSL mechanism with the adiabatic cooling of the same due to 
Hubble expansion. The slow cooling rate due to Hubble expansion ensures that any heating 
mechanism during that period should also be very small, constraining the parameter of the 
model to a much smaller value (as small as 10 -8 s _1 ). While Adler's analysis considers 
bounds yielding due to maintaining thermal equilibrium of the IGM, the analysis done in 
this article studies the breakdown of thermal equilibrium between photons and matter in 
cosmic plasma before recombination due to CSL heating and thus constrains the energy in- 
take by the photons in the cosmic plasma before z > 1100. As the processes, like Compton 
scattering, are very efficient in redistributing the energy over the spectrum and thus ther- 
malizing any anomalous heat-injection, a large enough heat injected through CSL can be 
tolerated by the cosmic soup without generating large spectral distortions and thus yielding 



a weaker bound on the strength parameter than the one obtained from IGM heating [13]. 
On the other hand, while calculating the heat intake by the CMBR photons due to heating 
of protons by CSL mechanism over the age of the universe, the case which is similar to the 
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ones considered in this article, Adler too found a bound of Ao < 10~ 5 s -1 [13( similar to the 
ones obtained here. 

It is nevertheless very important to mention at this point that the above bounds have 
been obtained by considering only CSL effect as the sole energy injecting process before 
recombination. However there can be other mechanisms which can lead to injection of en- 
ergies at such high-redshifts like adiabatic cooling of ordinary matter 3o[ 



black holes 



191 ]. decaying particles [l8|, |22], dissipation of magnetic fields 



31 1, evaporating 
32] or supercon- 



ducting strings [33] . Some of these processes are theoretically well understood and thus the 
distortion associated with such mechanisms can be obtained precisely. Hence the bounds on 
the strength parameter A obtained in this article are the most conservative ones. Any other 
energy-release mechanism can further tighten the bounds obtained here. Furthermore, in 
this article only white noise field has been considered for the stochastic field Wt(x) required 
for the collapse processes in CSL model. There have been attempts motivated from energy 
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361 ] . which is primarily 



conservation to model CSL mechanism with non-white noise 
dominant in low frequency regime. It could be interesting to analyze this kind of collapse 
mechanism in the cosmological context and can be ventured in future studies. 



Appendix A: Brief derivations of fi— type distortion 

Here we will give a brief derivation of fi— distortion following [ljlQ]. For a small chemical 
potential (fi <C 1) the energy density of a Bose-Einstein distribution is given as 

e = — — M-3— /i), (Al) 

where <tb ~ 5.67 x 10~ 8 kg s -3 K~ 4 is the Stefan-Boltzmann constant. The number density 
of photons in this distribution is 

n =^(l-2£,A (A2) 



kc V I 2 

The constants I n are related to the Reimann ( function as 

dx-^- [ = n\((n + l), (A3) 

with the numerical values I\ ~ 1.645, I2 — 2.404 and I3 ~ 6.494. Energy release processes in 
the early epochs may involve direct energy injection or injection of high frequency photons. 
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The change in the energy density given in Eq. (lAlj) is given by 



c j 

— = 4 51nT-3-^/x, (A4) 

£ ^3 



whereas the change in the number density given in Eq. (1A2|) will be 

— = 3 51nT-2-fy. (A5) 
n J 2 

Solving the above two equations simultaneously one obtains the change in chemical potential 
of the photons in the Bose-Einstein distribution as given in Eq. (j5T7|) . 



Appendix B: Brief derivations of y— type distortion 



Here we will give a brief derivation to calculate the y— parameter generated due to energy 
injection. In the non-relativistic limit, described by KbT e <C m e c 2 and k^T^ <C m e c 2 , Kom- 
mneets' equation is the kinetic equation of Compton scattering which can be expressed as 



22] 



<9n 7 
~dt 



N e a T c 



n b T e \ 1 d 



rripC 2 J xl dx f 



(Bl) 



where n 7 is the photon occupation number and x 



when T e is much larger, yielding 



hu 



The above equation simplifies 



N e a T c 



K, b T e \ 1 d 



x 



4 <9n 7 



(B2) 



dt c \m e c 2 J x\dx e \ e dx e 
Using the definition of Compton y— parameter given in Eq. (166]) the above equation can be 
written as 



1 d 



4 <9n 7 



(B3) 



dy x\ dx e \ e dx e 

where we have dropped the term T 7 in Eq. fltjol) as in this case T e > T T The energy density 
e 7 of photons can be written in terms of occupation number Ttry CIS 



87lh f 3 J 87F / 



Thus multiplying both sides of Eq. (IB 3D by x\ and solving for the integration yields 

fe 7 
5y 



4e 7 . 



(B4) 



(B5) 



Thus energy injection of fe 7 during 2<5x 10 4 can yield a Compton y— parameter as given 
in Eq. (jH 
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